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Evaluating multiloop Feynman integrals by Mellin-Barnes representation 

V.A. Smirnov''* 

'^ Nuclear Physics Institute of Moscow State University, 
Moscow 119992, Russia 

The status of analytical evaluation of double and triple box diagrams is characterized. The method of Mellin- 
Barnes representation as a tool to evaluate master integrals in these problems is advocated. New MB represen- 
tations for massive on-shell double boxes with general powers of propagators are presented. 
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1. Introduction 

When calculating physical quantities that de- 
scribe a given process one needs to evaluate a lot 
of Feynman integrals. After a tensor reduction 
based on some projectors (see, e.g., PP) a given 
Feynman graph generates various scalar Feynman 
integrals that have the same structure of the in- 
tegrand with various distributions of powers of 
propagators. A straightforward analytical strat- 
egy is to evaluate, by some methods, every scalar 
Feynman integral generated by the given graph. 
If the number of these integrals is small such 
strategy is reasonable. In non-trivial situations, 
where the number of different integrals can be at 
the level of hundreds and thousands, it is reason- 
able to follow a well-known advanced strategy: to 
derive, without calculation, and then apply inte- 
gration by parts (IBP) |2| and Lorentz-invariance 
(LI) PI identities between the given family of 
Feynman integrals as recurrence relations. The 
goal of this procedure is to express a general inte- 
gral from the given family as a linear combination 
of some basic (master) integrals. Therefore the 
whole problem of evaluation is decomposed into 
two parts: solution of the reduction procedure 
and evaluation of the master Feynman integrals. 

There were several recent attempts to make the 
reduction procedure systematic: 

(i) Using the fact that the total number of IBP 
and LI equations grows faster than the number 
of independent Feynman integrals one can sooner 
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or later obtain an overconstrained system of equa- 
tions msi. 

(ii) Using relations that can be obtained by 
tricks with shifting dimension [Sj. 

{in) Baikov's method [7]. 

(iv) Another attempt in this direction is the use 
of Grobncr basis (see, e.g., |H]). 

To evaluate master integrals one uses various 
methods. In particular, for the evaluation of dou- 
ble and triple boxes the following two methods 
were successfully applied last years: the method 
of Mellin-Barnes (MB) representation j5] and dif- 
ferential equations ^U] • The first of them is based 
on the following representation 
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which is applied to replace a sum of two terms 
raised to some power by their products in some 
powers at the cost of introducing an extra integra- 
tion. Usually one starts from alpha or Feynman 
parameters, then introduces, in a minimal way, 
MB integrations, performs internal integrations 
over Feynman parameters in terms of gamma 
functions and obtains a multiple MB represen- 
tation. It is useful to derive such representation 
for general powers of the propagators. 

The standard procedure of taking residues and 
shifting contours is used [H], with the goal to ob- 
tain a sum of integrals where one may expand in- 
tegrands in Laurent series in e = (4 — (i)/2 (where 
d is the space-time dimension within dimensional 
regularization }lll). Then one can use the first 



and the second Barnes lemmas and their corol- 
laries to perform some of the MB integrations 
explicitly. In the last integrations which usually 
carry dependence on the external variables, one 
closes contour in the complex plane and sums up 
corresponding series. (See ^] for details of the 
method) . 

A typical example of application of this method 
is the evaluation of master integrals for massless 
on-shell (pf = 0, i = 1, 2, 3, 4) double boxes pTS] 
where the reduction procedure was performed us- 
ing shifting dimension |13ll4j . with multiple sub- 
sequent applications |16| . 

2. Triple boxes and double boxes with one 
leg ofF-shell 

The first calculation of massless on-shell triple 
boxes in dimensional regularization was done in 
|17| . where the Regge asymptotics (in the limit 
t/ s — > 0) of the master planar triple box shown 
in Fig.|21was calculated with the help of the strat- 
egy of expansion by regions ^21 ■ Later it became 
possible to evaluate it analytically jJUj, with the 
help of a sevenfold MB representation, in terms of 
harmonic poly logarithms 20 . In fact, any mass- 




The complexity of the problem depends not 
only on the number of loops but also on the num- 
ber of parameters. One way to include one more 
parameter into massless on-shell double boxes is 
to consider one leg off shell, i.e. p\ = q^ ^ Q^ 
pf = 0, i = 2,3,4. This problem was successfully 
solved last years. The reduction to master inte- 
grals was done using Laporta's idea in |0 . Master 
integrals were calculated using MB representation 
(first results in TT) and DE (systematic evalua- 
tion in [5]). All results are expressed in terms 
of two-dimensional harmonic polylogarithms [H] 
which generalize harmonic polylogarithms. This 
combination of reduction based on jkb^ and DE 
was successfully applied in numerous calculations, 
e.g., various classes of vertex diagrams I23| . 

3. Massive on-shell double boxes 

Another way to include one more parameter 
into massless on-shell double boxes is to turn 
to massive on-shell double boxes with one non- 
zero mass, pf — m^, i — 1,2,3,4. Keeping in 
mind Bhabha scattering let us distinguish two 
planar and one non-planar diagrams as three ba- 
sic most complicated types of massive on-shell 
double boxes. The first planar graph is shown in 
Fig. [5] For the corresponding general Feynman 




Figure 1. Planar triple box diagram. 



less planar on-shell triple box can be evaluated 
by this procedure so that evaluation of three-loop 
virtual corrections to various scattering processes 
is not a miracle. 

Let us mention that the necessity to analyti- 
cally evaluate massless on-shell double and triple 
boxes arose when studying cross order relations 
in A'^ = 4 supersymmetric gauge theories |21|. In 
fact, to check some factorization one needs just 
one more triple box, in addition to Fig. 1. 



Figure 2. Planar massive on-shcU double box of 
the first type 
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with an irreducible numerator chosen as (fc+pi + 
P2+P3)^, the following sixfold MB representation 
was derived in |24| : 
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where 04557 = 04 + 05 + 06 + 07, 0122388 = ai + 
202 + 03 + 2o8, Z1123 = 2zi + Z2 + Z3, etc. 

Analytical evaluation of the master double box 
BpL,i{^, ...,1,0) was performed in P^, with a 
result in terms of polylogarithms. In [53], pla- 
nar 2-loop box diagrams with one-loop insertion 
were also by DE. The general MB representation 
(Pj) can be used for the evaluation of other mas- 
ter integrals. As an example, a double box with 
a numerator, BpLi{l, . . . , 1, —1) is evaluated in 
|26|. The finite part of the result in e includes 
polylogarithms and HPL 
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The second planar graph is shown in Fig. 13 
For the corresponding general Feynman integral 




Figure 3. Planar massive on-shell double box of 
the second type 
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with the same irreducible numerator, the follow- 
ing sixfold MB representation can be derived: 
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This multiple MB representation was used in 
|26| for the evaluation of the master planar dou- 
ble box BpL_2(l, • ■ • , 1,0; s,i,r7i^;e). The result 
includes polylogarithms, HPL as well as two- 
parametric integrals of elementary functions. It 
is not clear at the moment whether these inte- 
grals can be written in terms of HPL or 2dHPL 
depending on special combinations of s, t and m^, 
or a new class of functions is needed. This result 
as well as many other results for double and triple 
boxes mentioned in this short review is confirmed 
by numerical evaluation by means of a method 
based on a sector decomposition in the space of 
alpha parameters [77| . 

The non-planar graph is shown in Fig. ^ 
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Figure 4. Non-planar massive on-shell double box 



The following eightfold MB representation of 
the general non-planar double box, with the same 
irreducible numerator as above, can be derived: 
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This representation can be used for the evaluation 
of the non-planar master planar double boxes. 

4. Perspectives 

It has been reported |^ that the reduction 
of the massive on-shell double boxes relevant to 
Bhabha scattering can be done using Laporta's 
algorithm. Master integrals for these three re- 
cursion problems were identified. The calcula- 
tion of these master integrals was performed in 
all the cases when four indices o^ are positive 
and for some partial cases with five positive in- 
dices. In this problem the method of DE meets 
some complications because differential equations 
higher that the second order appear for compli- 
cated master integrals. Whether or not this ob- 
stacle can be overcome is an open question. How- 
ever, the three general MB representations pre- 
sented above can be certainly used for the ana- 
lytical evaluation of the master integrals. 

Let us finally characterize advantages of the 
method based on MB representation: 

(i) The MB representation can be derived for 
a general Feynman integral corresponding to a 
given graph, i.e. with general integer powers of 
the propagators. 

(ii) Resolution of singularities in e is much sim- 
pler than in alpha and Feynman parametric inte- 
grals. 

(m) After the resolution of singularities in e 
one can always switch to numerical evaluation, 
at least in order to check analytical results. The 
convergence along imaginary axis is always per- 
fect. 

[iv) Automation of calculations based on MB 



representation looks promising. 
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